APPLIED  PHYSICS  LETTERS  95,  031912  (2009) 


Extension  of  the  diffuse  mismatch  model  for  thermal  boundary  conductance 
between  isotropic  and  anisotropic  materials 
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This  model  is  an  extension  of  the  diffuse  mismatch  model  (DMM),  tailored  to  accurately  predict 
thermal  boundary  conductance  (hBD)  at  interfaces  where  one  material  comprising  the  interface  is 
characterized  by  high  elastic  anisotropy.  Temperature-dependent  specific  heat  is  calculated  with  this 
vibrational  model  and  compared  to  published  values.  Modifications  to  the  DMM  that  incorporate  the 
vibrational  model  are  presented  with  predictions  of  hBD  at  a  metal-graphite  interface.  This  model 
slightly  underestimates  experimental  data,  as  expected,  as  the  large  acoustic  mismatch  between 
metals  and  graphite  suggests  inelastic  scattering,  something  the  DMM  does  not  take  into  account. 
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The  ability  to  accurately  predict  thermal  boundary  con¬ 
ductance  (/?bd)  at  interfaces  is  essential  for  the  effective  de¬ 
sign  of  nanostructured  devices,  especially  thermoelectrics,1 
and  thermal  interface  materials ?  One  of  the  most  widely 
implemented  models  for  predicting  /?BD  is  the  diffuse  mis¬ 
match  model  (DMM).’  The  DMM  assumes  that  all  incident 
phonons  scatter  diffusely  and  elastically  at  the  interface  be¬ 
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tween  two  materials.  The  probability  of  a  phonon  diffusely 
scattering  across  the  interface  is  thus  related  only  to  the  oc¬ 
cupied  density  of  states,  the  phonon  group  velocity,  and  the 
phonon  energy  on  each  side  of  the  interface.  The  transmis¬ 
sion  coefficient  from  side  1  to  side  2,  which  describes  this 
probability,  is  given  by 


_ ^jJ^DJhwv2JD(at,  v2J)f  BEdco _ 

^3jfoDih(ovijD((o,  v\  j)/Be<7ui  +  Sjf^D4hQ)v2jD(to,  v2j)fw.dw  ’ 


where  j  is  a  particular  polarization,  at  is  the  phonon  fre¬ 
quency,  coD  is  the  cutoff  frequency,  v  is  the  phonon  group 
velocity,  D  is  the  density  of  states,  and  /Be  is  the  Bose- 
Einstein  distribution.  The  elastic-only  assumption  included 
in  the  development  of  the  DMM  enforces  that  the  upper  limit 
for  all  integrals  in  Eq.  (1)  is  the  cutoff  frequency  of  the 
material  with  the  lower  maximum  vibrational  frequency. 

Simplifications  to  £  can  be  made  for  interfaces  between 
isotropic  materials  such  that  the  expression  appears  in  a  more 
recognizable  form.  In  the  case  where  both  materials  that 
comprise  the  interface  are  three-dimensional  isotropic  Debye 
solids,  where  the  three-dimensional  Debye  density  of  states 
is  given  by 


D(a>,  v) 
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integration  over  at  in  Eq.  (1)  reduces  if  to 
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Integration  over  at  yields  a  frequency  independent  transmis¬ 
sion  coefficient.  Still,  the  form  of  t,  as  given  by  Eq.  (3)  is 
identical  to  the  form  presented  by  Swartz  and  Pohl.  Regard¬ 
less  of  the  formation  of  If,  the  DMM  yields  the  following 
expression  for  hBD : 


1  "  f"DJ  d 

h)fD1=-^2t  J  'C  '2h^V\,jD(ot,vUj)-^JBEdot.  (4) 

While  the  DMM  does  provide  an  estimate  of  hHl)  due  to  the 
difference  in  vibrational  properties  between  the  two  materi¬ 
als  that  comprise  an  interface,  it  can  largely  over  or  under¬ 
predict  experimentally  determined  values  for  many 
interfaces.43  Modifications  to  the  DMM  have  recently  been 
developed  to  account  for  both  mixing  at  interfaces,6  as  well 
as  the  contribution  of  inelastic  scattering,7  and  both  have 
shown  a  great  improvement  in  overall  accuracy  of  the  model. 
This  particular  extension  of  the  DMM  accounts  for  interfaces 
between  isotropic  and  highly  anisotropic  materials  at  el¬ 
evated  temperatures.  Example  calculations  are  given  for 
metal-graphite  interfaces. 
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FIG.  1.  (Color  online)  The  specific  heat  of  graphite  as  a  function  of  tem¬ 
perature  from  100  to  2000  K.  The  solid  line  is  the  value  calculated  using  the 
effective  density  of  states,  Z)eff,  and  the  data  points  are  published  values 
(Ref.  8).  As  expected,  this  model  underestimates  the  specific  heat  at  lower 
temperatures  as  a  result  of  ignoring  interlayer  vibrational  modes.  Still,  it  is 
evident  that  the  model  converges  at  the  high-temperature  limit. 


FIG.  2.  (Color  online)  Predicted  values  of  /;BD  as  a  function  of  temperature 
from  50  to  500  K  for  the  Al-graphite  interface  for  both  a-  and  c-axis  graph¬ 
ite  orientation.  The  predicted  value  for  the  a-axis  orientation  is  over  an  order 
of  magnitude  higher  than  for  the  c-axis  orientation  for  the  entire  temperature 
range. 


Graphite  is  a  highly  anisotropic  material,  with  a-axis 
(within  the  basal  plane)  and  c-axis  (perpendicular  to  the 
basal  plane)  thermal  conductivities  differing  by  nearly  three 
orders  of  magnitude.*  Any  appropriate  extension  of  the 
DMM  must  take  this  anisotropy  into  account.  While  other 
models  have  calculated  /zBD  at  metal-graphite  interfaces,9 
these  calculations  involve  derivation  and  integration  of  the 
complex  dispersion  relations  of  graphite.  However,  for  pre¬ 
dictions  of  hBD  at  metal-graphite  interfaces  at  temperatures 
above  cryogenic,  the  strong  anisotropy  of  graphite  allows  for 
certain  simplifying  assumptions.  Interlayer  bonding  in  graph¬ 
ite  is  very  weak  and  is  often  described  as  a  van  der  Waals 
type  interaction.  The  dispersion  diagrams  of  graphite  show 
that  the  acoustic  interlayer  vibrations  that  participate  in  ther¬ 
mal  transport  exist  at  frequencies  below  3  THz,  one  to  two 
orders  of  magnitude  below  the  maximum  longitudinal  and 
transverse  acoustic  vibrational  frequencies  within  the  basal 
plane.111  Additionally,  studies  looking  at  both  graphite  and 
carbon  nanotubes  have  indicated  that  these  interlayer  vibra¬ 
tional  modes  saturate  at  temperatures  above  50  K.  ' 

Previous  work  has  shown  that  as  a  result  of  this  weak 
interlayer  vibrational  coupling,  thermal  conductivity  of 
graphite  can  be  determined  by  treating  vibrations  as  a  one-  or 
two-dimensional  phonon  gas  depending  on  orientation.1"  A 
similar  reduction  can  be  made  when  calculating  other  ther¬ 
mal  properties.  At  elevated  temperatures,  the  number  of 
populated  vibrational  modes  within  the  basal  plane  far  ex¬ 
ceeds  the  number  of  interlayer  modes,  as  evident  by  the  dis¬ 
persion  diagrams  of  graphite,1013  and  the  low  temperature  at 
which  interlayer  modes  saturate.911  Thus,  the  vibrational 
spectrum  of  each  layer  can  be  described  by  a  two- 
dimensional  Debye  density  of  states, 


co 

D2d(u,v)  =  - — 9,  (5) 
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where  co  is  the  frequency  of  a  given  phonon.  To  determine 
the  density  of  states  of  an  assembly  of  monolayers,  this  two- 
dimensional  density  of  states  must  be  scaled  by  N,  or  the 


number  of  monolayers  per  unit  length.  Thus,  the  effective 
density  of  states  for  the  graphite  bulk  is  given  by 


Deff(w,  v)  =  N 


CO  CO  1 

27 tv1  2ttv^  d ' 


(6) 


where  d=  3.35  A  is  the  interlayer  spacing  for  graphite.14  Un¬ 
der  this  model,  the  single  monolayer  and  the  graphite  bulk 
have  the  same  maximum  vibrational  frequency;  the  only 
thing  that  differs  is  the  number  of  available  modes  at  each 
frequency.  Note  that  three  polarizations  are  considered  be¬ 
cause  the  assembly  of  two-dimensional  subsystems  creates  a 
three-dimensional  bulk  system.  While  the  vibrational  modes 
of  each  monolayer  or  subsystem  are  not  coupled,  forces  be¬ 
tween  the  subsystems  are  responsible  for  the  third  polariza¬ 
tion.  The  electronic  contribution  to  thermal  properties  is  not 
considered  in  the  following  derivation,  as  the  number  of  ef¬ 
fectively  free  electrons  within  the  basal  plane  has  been  re¬ 
ported  to  be  very  low  as  compared  to  that  of  a  typical  metal 
(neff=2.3  X  10-4  electrons  per  atom).13 

In  order  to  validate  this  treatment  of  the  vibrational 
modes,  the  specific  heat  of  graphite  was  considered.  Calcu¬ 
lations  of  the  specific  heat  of  crystals  with  a  basis  greater 
than  one  must  account  for  optical  phonons  via  the  Einstein 
model  in  addition  to  the  acoustic  phonon  contribution  using 
the  Debye  model.  The  formulation  of  the  specific  heat  of 
Debye  and  Einstein  solids  is  readily  available  in  other 
references.16  From  examination  of  the  vibrational  spectrum 
and  the  dispersion  diagrams  for  graphite,10’13  the  optical 
peaks  occur  at  the  acoustic  cutoff  frequencies  for  each  polar¬ 
ization.  The  cutoff  frequencies  for  each  polarization  are  de¬ 
termined  by 

ix;,4t7.V.  (7) 

where  N  is  the  planar  lattice-point  density  of  a  single  mono- 
layer.  Using  this  cutoff  frequency,  the  specific  heat  is  calcu¬ 
lated  from  100  to  2000  K. 

Figure  1  shows  a  comparison  of  specific  heat  calculated 
with  the  density  of  states  given  by  Eq.  (6),  the  specific  heat 
calculated  with  the  isotropic  Debye  density  of  states  given  by 
Eq.  (2)  using  both  a-  and  c-axis  properties  and  published 
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values.  The  specific  heat  calculated  using  Eq.  (6)  matches  the 
published  data  well  at  higher  temperatures  and  converges  to 
the  high-temperature  limit.  The  small  discrepancy  between 
calculated  and  reported  values  may  lie  within  the  assumption 
of  the  Einstein  model  that  all  optical  phonons  vibrate  at  the 
same  frequency.  Examination  of  the  dispersion  relationships 
reveals  that  this  is  not  entirely  accurate,  but  still  provides  a 
good  estimate.10  Since  this  model  is  concerned  specifically 
with  transport,  a  more  accurate  treatment  of  the  optical 
phonons  is  not  of  primary  interest  since  optical  phonons  do 
not  greatly  contribute  to  thermal  transport.  The  calculated 
values  of  specific  heat  using  the  three-dimensional  isotropic 
Debye  density  of  states  do  not  match  the  published  data, 
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demonstrating  that  an  isotropic  description  of  graphite  is  not 
accurate. 

When  using  the  effective  density  of  states  given  by  Eq. 
(6)  to  calculate  both  £  and  /zBD,  one  must  be  careful  to  con¬ 
sider  the  phonon  group  velocities  for  each  polarization  in  the 
direction  of  transport,  yet  only  use  the  in-plane  phonon  ve¬ 
locities  when  calculating  the  density  of  states.  Again,  if  each 
graphite  monolayer  is  treated  as  an  independent  two- 
dimensional  subsystem,  its  vibrational  states  are  dictated  by 
its  two-dimensional  properties,  or  under  the  Debye  approxi¬ 
mation,  the  in-plane  phonon  group  velocity.  Thus,  'C  '2, 
where  heat  flows  from  side  1,  a  metal  film,  to  side  2,  a 
graphite  substrate,  becomes 


_ Sy/oD-^Mi>2JDeff(M,yai)/BEr/fa) _ 

2//oDjfr Mv2jDeit( a,  vaj)fBEda)  +  2  Vij)fBEd(»’ 


at  which  point  hBD  can  be  calculated  by  Eq.  (4).  Figure  2 
shows  hBD  as  a  function  of  temperature  from  50  to  500  K  for 
the  Al-graphite  interface  for  two  different  orientations  of 
graphite.  The  predicted  value  for  the  a-axis  orientation  is 
over  an  order  of  magnitude  higher  than  for  the  c-axis  orien¬ 
tation  for  the  entire  temperature  range.  This  result  seems 
logical,  both  intuitively  and  conceptually,  as  the  thermal  con¬ 
ductivity,  and  hence,  the  phonon  group  velocity,  is  much 
higher  in-plane  compared  to  out-of-plane.  The  room  tem¬ 
perature  /?BD  at  the  Al-graphite  interface  (c-axis)  has  been 
reported  as  50  MW/m2  K.17  The  model-predicted  value  is 
30  MW/m2  K,  a  40%  underprediction.  This  is  to  be  ex¬ 
pected,  however,  as  the  elastic-only  assumption  made  by  the 
DMM  will  significantly  limit  the  number  of  phonons  that  can 
participate  in  transport  at  the  interface.  Still,  the  accuracy  of 
this  model  is  quite  good  when  considering  the  order-of- 
magnitude  error  often  associated  with  the  DMM.5 

Overall,  the  underlying  physics  of  the  graphite  structure, 
as  well  as  a  comparison  of  predicted  specific  heat  to  pub¬ 
lished  values,  establishes  the  validity  of  the  model  when  cal¬ 
culating  thermal  properties  at  elevated  temperatures.  An  ex¬ 
tension  of  the  DMM  employing  this  vibrational  model  has 
been  presented,  with  applicability  to  predicting  hBD  at  inter¬ 
faces  between  isotropic  and  highly  anisotropic  materials. 
Comparison  between  the  predicted  and  experimental  room 
temperature  hBD  at  an  Al-graphite  interface  shows  good 
agreement,  and  the  underprediction  is  to  be  expected.  This 
model  can  easily  be  extended  to  other  anisotropic  media 
where  weak  coupling  between  subsystems  is  present.  Further 
experimental  data  will  be  required  to  evaluate  the  tempera¬ 
ture  range  over  which  this  model  is  valid. 
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